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Abstract 

' The 'minimal' payment — a payment method which minimizes the number of coins in a 

. purse — is presented. We focus on a time series of change given back to a shopper repeating 

5— ( \ the minimal payment. The delay plot shows visually that the set of successive change 

possesses a fine structure similar to the Sierpinski gasket. We also estimate effectivity of 
the minimal-payment method by means of the average number of coins in a purse, and 
\ conclude that the minimal-payment strategy is the best to reduce the number of coins 

in a purse. Moreover, we compare our results to the rule-60 cellular automaton and the 
Pascal-Sierpinski gaskets, which are known as generators of the discrete Sierpinski gasket. 
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An English proverb says that "a heavy purse makes a light heart" , but we do not welcome a 
J> ■ literal heavy purse full of coins. By repeated payments in everyday shopping, the number of 
coins in the shopper's purse constantly increases or decreases, and sometimes the purse gets 
CN ! full of coins. A purse with a lot of coins is uneconomic because it becomes heavy and bulky, 
and even worse, the shopper has a lot of trouble in picking out coins from the purse. Electronic 
money [1] is hoped to eliminate such trouble, but it has not become a substitute for real cash 
so far because of technical, economical, and legal issues. 

Some mathematical problems treat of coins. A famous one is known as the change-making 
problem [2]: how can a given amount of money be made with the least number of coins of given 
^ \ denominations? This problem is a variation of the knapsack problem [3], and it is known to be 
NP-hard [1]. The greedy algorithm [5] gives a solution in some limited cases, and other heuristic 
algorithms have been also proposed [0H5]. The optimality or efficiency of a monetary system 
itself has been also evaluated P^TT] , and the Frobenius coin problem asks what is the 

largest amount of money that cannot be obtained using only coins of specified denominations. 
We note that such problems are for retailers (how to make change) or for governments (what 
currency system is preferable); consumers are not considered. To prevent the swelling of a 
purse, which is a problem for consumers, we propose a method to minimize the number of coins 
in the purse through the payment. We also show that this method accompanies a remarkable 
mathematical structure like fractals. 
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2 What is the minimal payment? 



We define the 'minimal payment' as follows, (i) A retailer keeps all kinds of coins sufficiently, 
and gives back change with the least number of coins. Namely, the retailer never gives 10-cent 
change by ten 1-cent coins or two 5-cent coins, (ii) The shopper pays with a combination of 
coins that minimizes the number of coins in the purse after the payment, (iii) The coins paid 
by the shopper and the ones given back to the shopper do not include common coins. 

Let us play with a simple example. Imagine a nation where the coins are issued in 06, <}2, 
and 01 denominations ('0' is the currency sign of this nation, and we also use this symbol as 
the generic currency sign in this paper, instead of the formal symbol 'D'), and a shopper in this 
nation has each of 06, 02, and 01 coins in the purse, and is buying an item of 05. We compare 
two ways of payment. (1) If the shopper pays with one 06-coin alone, then the change is 01, 
hence three coins (one 02- and two Ol-coins) are left in the purse. (2) If the shopper pays with 
one 06- and one Ol-coins, then the change is 02 (paid with one 02-coin), hence two 02-coins 
are left after the payment. Other payment patterns do not meet the condition (iii). Therefore, 
the minimal payment in this case is that the shopper pays with one 06- and one Ol-coins and 
receives one 02-coin in change (see Fig. [T]). The important fact is that the number of coins left 
in the shopper's purse can vary according to the combination of coins in the payment. 

Next we present the efficient calculation method of the minimal payment. Let us consider 
again the currency '0'- The notation (ni, 77.2, tt-s) represents the content of the shopper's purse, 
where rii, n2, and 77.3 are respectively the number of 06, 02, and 01 coins. In the above 
example, the coins before the payment are written as (1, 1, 1) (see also Fig. [T]). By using this 
vector-like notation, the total amount of money in the purse can be expressed as the scalar 
product: (1, 1, 1) ■ (06, 02, 01) = 1 ■ 06 + 1 ■ 02 + 1 ■ 01 = 09. In this sense, we identify a vector 
with the total amount, and simply express as 09 = (1, 1, 1). Returning to the above example, 
04(= 09 — 05) is left in the purse after the purchase of an item of 05. Through the transition 
from 09 = (1, 1, 1) (before the payment) to 04 = (0, 2, 0) (after the payment), one 06- and one 
Ol-coins are gone and one 02-coin is gained, which correspons to the minimal payment process 
(Fig. [1]). We note that this calculation procedure for the minimal payment is far more efficient 
than the 'full search', and it can be applied straightforwardly to general currencies. 

3 Fractality of successive minimal-payment processes 

Here we show the mathematical structure of successive minimal-payment processes. We set 
the following four assumptions, (i) The shopper and the retailer iterate the minimal payment, 
(ii) Prices are generated randomly and independently in each payment, and we used a uniform 
random number in the numerical calculations, (iii) In order for the shopper not to lack money, 
we assume that the shopper has sufficiently large number of banknotes, and that the shopper 
can use banknotes if needed, (iv) Every coin is a multiple of the immediately smaller coin, 
and the smallest banknote is a multiple of the largest coin. We call (iv) the 'multiplicable 
condition.' For example, the Japanese yen (1000 yen bill; 500, 100, 50, 10, 5, 1 yen coins), the 
Korean won (1000 won bill; 500, 100, 50, 10, 5, 1 won coins), and the Swedish krona (20 krona 
bill; 10, 5, 1 krona, 50 ore (= 0.5 krona) coins) are multiplicable. However, the majority of 
currencies in the world is partially not multiplicable; in the US dollar (1 dollar (= 100 cent) 
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bill; 50, 25, 10, 5, and 1 cent coins), 25 is not a multiple of 10. 

For the analysis of successive change in the minimal payment, we employ the delay plot. 
For a given time series z{0), z{l), ■ ■ ■ ,z{T), the delay plot is defined as a planar point set 
{{z{t), z(t + l))\t = 0, 1, ■ ■ ■ ,T — 1}. The delay plot is a useful tool for visualizing the correlation 
of adjoining pairs of the time series. In chaos theory [15] and nonlinear time series analysis [16], 
the delay plot is used for the reconstruction of an attractor from a series of experimental data. 
In general, a delay plot generated from a deterministic system forms an organized structure, 
whereas a plot from a stochastic noisy system forms an indistinct distribution. 

Figure [2] illustrates the delay plot of a series of change generated by computer simulation of 
the minimal payment. In the simulation, we adopted the Japanese yen as the currency system 
(which holds the multiplicable condition), and iterated payment processes 10'' + 1 times. The 
whole structure is similar to the celebrated fractal object 'Sierpinski gasket' [T7] . 

In order to explain the mechanism of this quasi-fractal structure, let us go back once again 
to the model currency '0', where 06, 02, and 01 coins are issued. We additionally imagine 
that the smallest banknote is 012 denomination. Note that this currency is multiplicable. 
(riiit) , n2{t) , n^ij:)) represents the coins in the shopper's purse after the t-th payment, and we set 
the initial condition (?t,i(0), n2(0), 77.3(0)) = (0,0,0). We also denote by c(t) = {ci{t) , C2{t) , c^it)) 
the change given back in the t-th payment. Since the prices are randomly generated, (ni, n2,ns) 
and (ci, C2, C3) are random variables. If ni{t) > ni(t — 1) (06 coins in the purse increases), the 
t-th change includes ni(t) —ni{t — l) pieces of 06 coin. Otherwise, if ni(t) < ni(t — 1), the t-th 
change includes no 06 coins. Hence we can simply express as ci(t) = max{ni(t) — ni(t — 1), 0}. 
The sum of the number of 06 coins in t-th and {t + l)-th change is then calculated as 

ci(t) + ci(t + 1) = max{ni(t) - ni{t - 1), 0} + max{ni(t + 1) - ni{t), 0} 

= max{ni(t + 1) - ni{t - l),ni{t) - ni{t - l),ni(t + 1) - ni(t),0}, (1) 

where we use the formula in a max-plus algebra [T^j: max{a, 6} + max{c, li} = max{a + 
c,a + d,b + c,b + d}. By the definition of the minimal payment and above algorithm, we 
can derive < ni(t) < 1 for any t > 0, where '1' at the right hand side is given by the 
ratio of two denominations as 012/06 — 1. Therefore, the four arguments of 'max' in Eq. 
([1]) are all equal to or less than 1 constantly, and we obtain < ci{t) + ci{t + 1) < 1 as a 
result. Similarly, we can derive the relations < C2(t) -|- C2(t + 1) < 2(= 06/02 — 1) and 
< C3(t) + C3(t + 1) < 1(= 02/01 — 1) for C2 and C3 respectively. These relations signify that 
successive change c(t) and c(t + 1) are correlated. 

By the above calculation, the points (c(t), c(t + 1)) of the delay plot are confined within the 
lattice points satisfying the conditions 

<ci(t) + ci(t + l) < 1, 

0<C2(t)+C2(t + l) <2, (2) 
< C3(t) + C3(t + 1) < 1. 

As in Fig. [3l the effect of including smaller coins progressively generates a nested pattern 
(smaller coins contribute to higher resolution). In each iterative step, the blocks (the squares 
with bold contours) are replaced with smaller blocks aligned in a step-shape (or an L-shape), 
and finally the structure of the delay plot has a fine structure as in the Sierpinski gasket. Note 
that the side lengths, 12 (the length of the 'whole' block), 6, 2, and 1, of blocks in the iterative 
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steps are identical with the denominations 012, 06, 02, and 01- This scenario of 'successive 
replacement' universally works on general multiplicable currency systems. In particular, if the 
denominations of coins grow geometrically (i.e., 01, 0^, 0^^, ■■■)) t^ie delay plot possesses 
(discrete) self- similarity, and its fractal dimension D is given by 

In ±±11 In 
D = ^-^ = l + '^. (3) 
m r m r 

This expression also appears in a study of Pascal's triangle [18]; we discuss delay plots and 
Pascal's triangle in more detail in Section |6l The multiplicable condition is necessary for the 
realization of a regular shape of the delay plot. For example, the delay plot of the US dollar 
exhibits a relatively disordered structure, reflecting the breaking of the multiplicable condition 
between 25 and 10 cent coins (Fig. H]). 



4 EfFectivity of the minimal payment 

Here we comment on the effectivity of the minimal-payment method, measured by the average 
number of coins in the shopper's purse. We estimated this average numerically. In this calcu- 
lation, we adopted the US dollar as the currency, and iterated the minimal payment 10'' times. 
The results are summarized in Fig. [51 Figure |5] (a) shows the fluctuation of the number of 
coins in the purse of the first 1000 steps, and Fig. [S] (b) is the numerically-generated probability 
distribution of the number of coins in the purse. By taking the mean value of Fig. (b), the 
number of coins in the purse is 4.20 pieces on average. Furthermore, the maximum number of 
coins in Fig. [5] (b) is eight. Thus, by using the minimal-payment strategy, we can constantly 
keep the coins in a purse 8 pieces or less. 

We can also calculate this average more rigorously. Since the prices are distributed randomly 
and independently in our numerical model, a probability of the amount of money in the purse 
becomes uniform. Figure |6] illustrates how many coins are needed to represent each amount of 
money, and the average number of coins is exactly 4.2 pieces by taking the mean value. Note 
that the two distributions in Fig. [5] (b) and Fig. [6] (b) are so alike. 

The average number can be calculated easily for a multiplicable currency. For the model 
currency '0' (012 bill; 06, 02, and 01 coins), the set of possible amounts of money in a purse 
is 

00 = (0,0,0), 01 = (0,0,1), 02 = (0,1,0), 03 = (0,1,1), 

04 = (0,2,0), 05 = (0,2,1), 06 = (1,0,0), 07= (1,0,1), (4) 

08 = (1,1,0), 09 = (1,1,1), 010 = (1,2,0), Oil = (1,2,1). 

For the same reason as above, each amount of money listed in Eq. emerges with a uniform 
probability. In Eq. (jl]), the possible numbers of 06 coins are either or 1, which appear six 
times each. Thus, the average number of 06 coins is 0.5, which is also given by (012 /06 — 1)/2. 
The same estimation is valid for the numbers of 02 and 01 coins. Therefore, the average 
number of coins in a purse is 

012/06- 1 06/02- 1 02/01 - 1 _ 2 
2 ^ 2 ^ 2 ~ ■ 
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We can straightforwardly apply this formula to any multiplicable currencies. For example, the 
average number of the Japanese yen (1000 yen bill; 500, 100, 50, 10, 5, and 1 yen coins) is both 
7.5 pieces. 

Conversely, we note that this average number can be regarded as the lower bound, that is, 
we cannot keep down the average number of coins below this lower bound, no matter what 
payment strategy is used. This is the 'minimality' of the presented payment method. 

5 Relevance to the rule-60 cellular automaton 

In the present section, let us consider the denominations of coins grow with a binary scale, that 
is, the coins are 01, 02, 04, 08, ■ ■ ■ , (}2"'~^, and the smallest banknote is 02". As in Section 
El let Ci{t) denote the number of 02"~* coins in t-th change, and the condition Cj(t) = or 1 
holds for any t and i. The corresponding delay plot is obtained by solving 

' < ci(t) + ci(t + 1) < 1, 

<C2(t)+C2(t+l) < 1, 

. (5) 

^ < Cnit) + Cnit + 1) < 1. 

In other words, c„(t) = 1 and c„(t + 1) = 1 never take place at the same time. The delay plot 
forms the discrete Sierpinski gasket (Fig. [7]). 

The discrete Sierpinski gasket appears in various systems in mathematics and physics, rang- 
ing from an analysis of the Tower of Hanoi [20] to pulse dynamics in a reaction-diffusion sys- 
tem [21] and a model of fracture [22]. As a typical example, we focus on the rule-60 cellular 
automaton in this section. 

A cellular automaton (CA) [23] is a discrete dynamical system on cells aligned linearly. The 
possible state of each cell is either '0' or '1', and the evolution of states generally depends on 
their own state and the states of their neighbors. In particular, the rule-60 CA can generate 
the discrete Sierpinski gasket, which is driven by the following rule, (i) If a cell and its left 
neighbor are in the same state, the next state of the cell is 0. (ii) Otherwise (i.e., the states 
of the two cells are different), the next state is 1. This evolution rule is schematically shown 
in Fig. [H] (a). We denote by the state of n-th cell at time t. Then, the evolution of the 
state is simply expressed as s^^^ = © -^^-i, where © is an exclusive-or operation defined as 
0©0 = 1©1 = and 0©1 = 1©0 = 1. A resultant pattern is depicted in Fig. M (b), 
where the initial condition is that one cell is '1' and the others are '0'. The evolution rules of 
other cellular automata, such as rule-90, 22, 154 and 102, can be expressed by ©, and they also 
generate discrete Sierpinski gaskets. 

Similarity between a delay plot of binary coins (Fig. [7]) and the evolution pattern of the 
rule-60 CA (Fig. [8] (b)) is obvious. We discuss this similarity more in detail. Figure [9] shows 
iterative steps for solving Eq. ([5]). As a condition for a smaller coin is considered, the blocks 
in each step are subdivided according to the following rule: 

(6) 

As shown in Fig. [10] (a), this subdivision preserves the exclusive-or structure of the rule-60 
CA. Moreover, a spatio-temporal pattern after the subdivision is similar to the one before 
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the subdivision, that is, the subdivision works on the rule-60 CA as a spatio-temporal scale 
transformation. Both the unit lengths of space and time are halved by subdivision (see Fig. 
ITU] (b)). Therefore, we conclude that the key to similarity between the delay plot of a binary 
currency and the rule-60 CA is compatibility of the rule ([6]) of subdivision and the exclusive-or 
relation. 

6 Comparison with the Pascal-Sierpinski gaskets 

Pascal's triangle (Fig. [TT] (a)) also generates a discrete Sierpinski gasket by coloring the odd 
numbers (Fig. [TT] (b)) [21]. Clearly, Fig. [TT] (b) is essentially the same as Fig. [7] and Fig. [8] 
(b). In fact, the coloring rule is reduced to the exclusive-or operation [25]. As a generalization, 
the Pascal-Sierpinski gasket [2S] of modulo r is obtained by coloring the elements not divisible 
by r. In this section, we compare the Pascal-Sierpinski gasket of modulo r and the delay plot 
of the coins 01, 0?", 0^^, ■ ■ ■ . 

In Fig. [12] (a), we show the delay plots corresponding to r = 3,4,5, and 6. The shapes 
of them are quite similar to the Sierpinski gasket. In contrast, the Pascal-Sierpinski gaskets 
are more complicated (Fig. [I2](b)). In the case of r = 3 and 5, the Pascal-Sierpinski gaskets 
are similar to the corresponding delay plots. However, the Pascal-Sierpinski gaskets of modulo 
r = 4 and 6 exhibit disordered structure, and they looks different from the delay plots of such 
r. According to the theory [27] of the Pascal-Sierpinski gasket, the Pascal-Sierpinski gasket 
possess the same structure as the corresponding delay plot if r is a prime number; however, the 
pattern is more complicated if r is a composite number. 

Analogues of Fig. [H] Eq. (|3|), and Eqs. ^ and ^ have appeared in some research articles 
of Pascal's triangle [181 [28] [29] . Hence, we note that the minimal payment is a real- world 
counterpart of such purely mathematical objects. 

7 Discussion 

In the numerical simulations of the minimal payment in this paper, we assume that the prices 
are distributed randomly, independently, and uniformly. Randomness and independentness are 
reasonable for real payment, but uniformness is a debatable point. In fact, the range of prices 
depends on each person's life style or consumption behavior. We simply model this non-uniform 
occurrence of prices by using a triangular random number. In this numerical calculation, the 
Japanese yen (1000 yen bill; 500, 100, 50, 10, 5, and 1 yen coins) was used as a currency, and we 
iterated the minimal-payment processes 10 ''-f 1 times. The prices are generated by a triangular 
random number on the interval [0, 1000], and the position of the peak of the triangle is a control 
parameter (Fig. [T3] (a)). The result is shown in Fig. [T3] (b). The average numbers of coins in 
the purse are within 7.5 ± 0.003, which is irrespective of the position of the peak of triangle. 
Thus, the average number of coins in a purse is still 7.5 pieces also in this more realistic case. 
We expect that the average number of coins in a purse remains constant over a wide class of 
probability distributions of prices. 

The term 'fractal time series' or 'temporal fractal' conventionally means that a time series 
has scale-invariance or power-law spectrum [5U] - [32] . On the other hand, the time series of 
change in successive minimal payment is another type of 'fractal time series' in the sense that 
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its delay plot exhibits self-similarity. We note that the power spectrum of a time series of 
change does not exhibit power law. 

The vector-like notation (ni,n2,n3) is a generalization of the binary (or r-adic) expansion, 
which is an essential tool in computer science [33l[3l]. A cellular automaton was invented as a 
prototype of digital computing. Namely, our analysis originates from analogue payment, and 
leads to digital fractals. In the analysis, we assume that a retailer gives back change with the 
least number of coins. However, real retailers may not give change with the least number of 
coins because they may run out of coins of certain denominations, or may feel that finding the 
least number is troublesome. This causes the increase of the number of coins in a purse. Thus, 
realization of the 'minimal' requires the retailers' cooperation, but we emphasize that we can 
avoid immoderate swelling of a purse at least if we pay attention to our payment styles. 

8 Conclusion 

We demonstrated that a payment method affects the number of coins left in a purse, and showed 
how to minimize the number of coins in a purse after a payment. We called this method the 
'minimal payment'. By using the delay plot, a time series of change under successive minimal- 
payment processes a fine structure like the Sierpinski gasket. The mechanism of this fine 
structure is reduced to a coupled inequalities as in Eqs. ([2]) and ([5]). We also compared delay 
plots and the rule-60 cellular automaton and the Pascal-Sierpinski gasket. 
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not minimal minimal 



Figure 1: Illustration for the minimal payment in a simple situation: coins are ()Q, 02, and 
01. The notation (1, 1, 1) means that one 06-, one 02-, and one Ol-coins are present in the 
purse. A shopper is buying an item of 05. Compare two types of payment (left: not minimal, 
right: minimal). 
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Figure 2: The delay plot of a series of successive change in the Japanese yen (1000 yen bill; 500, 
100, 50, 10, 5, 1 yen coins), constructed from a minimal-payment simulation iterated 10^ + 1 
times. Inset: the Sierpinski gasket. 
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Figure 3: The iterative steps for the determination of the possible area of the delay plot in the 
case of the currency '(>' (012 bill; 06, 02, and 01 coins). The pattern gets higher resolution as 
conditions for smaller coins are included. The unit squares are used to indicate lattice points. 
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Figure 4: The delay plot of a non-multiplicable case: the US dollar (1 dollar bill; 50, 25, 10, 
5, and 1 cent coins). 
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The number of coins in tlie purse 
(a) (b) 



Figure 5: (a) The time siries of the number of coins in the shopper's purse, generated by a 
numerical simulation of the minimal payment. The currency is the US dollar (1 dollar bill; 50, 
25, 10, 5, and 1 cent coins), and the first 1000 steps are shown, (b) The numerically-computed 
probability distribution of the number of coins in the purse. The average of this distribution is 
4.20 pieces, indicated by the dashed line. 




Figure 6: (a) The least number of coins that the US dollar needs to represent each amount 
of money, (b) The distribution of the number of coins, built from (a). The average number of 
coins is 4.2 pieces, indicated by the dashed lines. 
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Figure 7: A delay plot of a binary currency (the coins are 01, 02, 04, 08; and the smallest 
banknote is 016) forms a discrete Sierpinski gasket. 



space 



I 

□ 



□ 







(a) 



(b) 



Figure 8: (a) Evolution of the rule-60 cellular automaton. White cells represent the state 
'0', and colored ones represent the state '1'. (b) A spatio-temporal pattern of the rule-60 
cellular automaton. Horizontal direction represents space, and vertical direction represents 
time. Initially, one cell is the state 1, and the others are 0. 
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Figure 9: The first three iterative steps for solving Eq. (|5]). We can reach the solution by 
simply applying the rule ([6]) to each block recursively. 
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Figure 10: (a) Subdivision of the elementary evolution rule of the rule-60 cellular automaton. 
The small cells again possess the exclusive-or structure, (b) Subdivision of a spatio-temporal 
pattern of Fig. |8] (b), acting as a scale transformation. 
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Figure 11: (a) A part of Pascal's triangle. The elements are aligned on the square lattice, (b) 
The discrete Sierpinski gasket is obtained by coloring the odd numbers. 
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Figure 12: (a) The delay plots of the coins growing geometrically (01, ()r^, ■■■) for 
r = 3,4,5, and 6. They all exhibit ordered shapes like the discrete Sierpinski gasket. Since 
we compare only the shapes, the scales are omitted from the plots, (b) The Pascal- Sierpinski 
gasket of modulo r for the case r = 3,4, 5, and 6. The upper part of 100 rows are shown. 
Pascal- Sierpinski gaskets are similar to the delay plots for r = 3 and 5, and are different for 
r = 4 and 6. 
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Figure 13: (a) The probability density function of a triangular random number. The currency 
is the Japanese yen, so that prices are distributed between yen and 1000 yen. The height 
is given by normalization, (b) A numerical result of the average number of coins in a purse, 
where the prices are generated by a triangular random number. The average numbers are 
within 7.5 ± 0.003, irrespective of the position of the peak of a triangular random number. 
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